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We study a two-dimensional system of bosons interacting with a fluctuating U(l) gauge field 
with overdamped dynamics. We find two instabilities of the condensed phase at T — 0: one to 
phase separation and another to a homogeneous non-superfluid (Bose metal) . The presence of both 
instabilities in the model is dependent on the low-energy form of the gauge field propagator. We 
discuss the relevance of our findings to the U(l) gauge theory of the t—J model. 
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I. INTRODUCTION 

A Bose liquid naturally becomes a superfluid at low 
temperatures. The existence of a non-superfluid liquid 
state (the "Bose metal") down to temperature has 
been a subject of recent controversyaBn. Such a scenario 
might ar i§9iiji-t|k c bosons interact with fluctuating exter- 
nal fieldsa'QErQu. A fluctuating vector potential couples 
directly to the gradient of the phase of the superfluid 
order parameter. This makes it a natural choice for de- 
stroying long-range phase coherence in a superfluid. 

In this paper, we study how a vector potential with 
overdamped dynamics might destroy superfluidity. Our 
choice xjf-overdamped dynamics is motivated by the gauge 
theoryod of the t-J model, relevant to the cuprate super- 
conductors. The existence of a Bose metal far below the 
degeneracy temperature is crucial for the picture of spin- 
charge separation in that theory. 

Our study is the first to consider the damping of the 
zero-point fluctuations of the gauge field. We find that 
a metallic phase may exist in the part of the phase di- 
agram relevant to the cuprates. Related models with 
different dynamics have been studied before. Using a 
quasistatic approximation for the gauge field, Lee, Lee 
and KirnL] showed the loss of superfluidity in a quantum 
Monte Carlo simulation. On the other hand, a metal- 
lic phase was not found with a propagating ( "MaxweU" ) 
gauge fieldEl. Instead, the system is unstable to phase 
separation (unless it is stabilised by long-ranged repul- 
sion.) 

In our treatment, we use the perturbative corrections 
calculated by Feigelman et ala in a simple renormalisa- 
tion scheme. In Ref. |, a self-consistent scheme was used 
to calculate the superfluid density in the presence of fluc- 
tuating magnetic fields. It was clear from that work that 
the superfluid response depends on the compressibility of 
the Bose liquid, and vice versa: the superfluid expels flux 
over a penetration depth, but that inevitably reduces the 
local density, exciting phonons and causing an apparent 
increase in the compressibility. Since there is the possibil- 
ity of phase separation (crudely speaking, the bosons and 
magnetic flux expelling each other), we want to track this 
inter-dependence in more detail. Our scheme resembles a 
frequency-dependent mean-field theory. We successively 
integrate out high-frequency fluctuations, generating at 



each stage corrections to the low-energy response func- 
tions — the superfluid density n s and the compressibility 
k. We believe that this will give a better picture of the 
instabilities of the superfluid phase. 

The form of the paper is as follows. We start by re- 
viewing the method by which charge dynamics of the t-J 
model can be obtained from a model of bosons interact- 
ing with a U(l) gauge field. We then calculate correc- 
tions to the compressibility, phase stiffness and Meissner 
response, in a perturbation expansion about the super- 
fluid state. We use these corrections in a renormalisa- 
tion scheme, which we apply to systems with both over- 
damped and propagating gauge fields. Finally, we discuss 
our results and make some conclusions about the effect 
of different gauge fields on the bosonic systems. 

II. THE MODEL 

We will now describe in more detail our model, in par- 
ticular, its motivation from the gauge theory of the t— 
J model. In this paper, we focus on two-dimensional 
bosonic systems with (imaginary-time) Lagrangian den- 
sities of the form: (Ti = c = e = 1) 
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where b x is a complex bosonic field with local repulsion 
U and mass m, and is a (two-vector) gauge field with 
propagator D. We label both spatial and temporal coor- 
dinates by x = (r,r). 

As already mentioned, we study the gauge fields that 
arise from the U(l) gauge theory of the t-J model. The 
t—J model describes electrons with hopping integral t and 
antiferromagnetic exchange J. The Coulomb repulsion 
is so strong that we assume that there is at most one 
electron per site. The Hamiltonian is: 
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with a constraint of no double occupancy at each site. 
The operator c\ a creates an electron at site i with spin 
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<7, n, is the number of electrons on a site, and S, = 
c ia< T a/3 c i/3 is the spin operator at a site. The sums 
are taken over pairs of nearest neighbours on a square 
lattice. nn 

A "slave boson" treatmentErtl of the constraint of sin- 
gle occupancy reduces this model to a Lagrangian den- 
sity of the form in equation ^. The electronic operator 
is split into a spin-half fcrmion and a charge-e boson: 
c ia = b\fi a . The single-occupancy constraint means the 
fermion and boson currents must be equal and opposite: 
Jf = — Jf,. The gauge theory enforces this constraint 
by making the boson and fermion currents interact indi- 
rectly via a vector potential. 

After this (exact) decoupling, we may obtain an effec- 
tive theory for the bosons. The derivation of this (ap- 
proximate) theory is shown in appendix [A], The effective 
Lagrangian density for the bosons takes the form of equa- 
tion [j]; the theory aims to capture the strong correlations 
in the t — J model, which arise from the constraint of no 
double occupancy. 

As explained in appendix |A|, the scalar part of the 
gauge field that arises in this treatment has been ab- 
sorbed into the local repulsion U. The boson density is 
the charge doping away from half-filling. This appears to 
be consistent with the charge carrier density from trans- 
port measurements. So, we will concentrate on boson 
densities small compared to integer filling of the lattice 
and we can ignore the possibility of a Bose Mott insulator 
in our system. 

The excitations of the fermionic subsystem can ex- 
change energy and momentum with the bosons. Integrat- 
ing out the fermions amounts to treating them as an effec- 
tive medium through which the gauge field propagates. 
In fact, the Lagrangian for the gauge field (last term in 
(|l|)) is determined by the dynamics of the particle- hole 
excitations of the fermions. It can be showru that the 
propagator, D UJntq , is directly related to the fermion cur- 
rent response functions. This is consistent with the con- 
straint on the boson and fermion currents. 

Let us choose the gauge V • a = so that only the 
transverse component, a -1 , is non-zero and D^ v has only 
one degree of freedom which we call D. As shown 
in appendix the spectral density of this propaga- 
tor is A(ui,q) = 2 Im[— l/ITg(w, g)] where ±lj£ = (S^ u — 
( l/j. < lv/q 2 )[i(Jf'J{}S,q ~ {nf/m{)5^} is the retarded trans- 
verse polarisation of the fermions. 

For simplicity, we will model the fermions as an 
isotropic Fermi liquid with mass rrif and a density rif 
which we choose to be close to the half-filling density 
on a lattice model. For the cuprates, this is expected 
to be reasonable near optimal doping. In other words, 
we ignore the "pseudogap" spin physics in underdoped 
materials. 

With this choice of isotropic fermions, the polarisation 
can be derived analytically. This saves a great deal of 
computational time in the numerical calculations carried 
out in the next section. However, it ignores the under- 
lying lattice. For this reason we will introduce a high- 



energy cutoff, fl, in the spectral density of the gauge field 
propagator. We choose ft = Ep so that particle- hole ex- 
citations with energies greater than the Fermi energy Ep 
are excluded — these would have involved quasiparticles 
outside the Brillouin zone of the underlying lattice. The 
form of the polarisation and the effect of the cutoff fl are 
discussed in appendix^. Taking all this into account, wc 
obtain: 
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At small frequencies and momenta this reduces to the 

form D Uniq = (a q \u; n \ + xq 2 ) where \ = (127rm f ) _1 
is the diamagnetic susceptibility of a two-dimensional 
Fermi liquid and a q — 2ri{/kpq is the conductivity due 
to Landau damping. (However, we do not work with this 
simple form because we find that strong scattering occurs 
at large wavevectors.) 

Equation H completes our definition of the model in (|l|) 
by specifying the spectrum of gauge field fluctuations. 
In the following sections we discuss instabilities of the 
superfluid state in this model. 



III. THE RENORMALISATION SCHEME 

In this section, we investigate instabilities of the super- 
fluid state in the general model given by equation |l|. We 
discuss the dependence of these instabilities on the form 
of the gauge field propagator. The interaction between 
the bosonic excitations and the gauge field affects the 
compressibility of the bosons and their phase stiffness, as 
well as the Meissner response of the gauge field. These 
corrections can be calculated in perturbation theory, but 
near instabilities of the superfluid state a more advanced 
treatment is necessary, so we introduce a simple renor- 
malisation scheme. 

We first expand about a uniform condensate by writ- 
ing: \P x = yjn + p x e l8j: , where n is the mean boson den- 
sity, p represents density fluctuations, and 6 represents 
fluctuations in the phase of the bosonic field. The La- 
grangian density becomes C x = C x + C" 
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The propagator (D°_ x ,)^ = (%6(x - x>) + D^,) 
where the first term is the Meissner response which mod- 
ifies the bare propagator. 

The Lagrangian density C° describes a superfluid Bose 
liquid with compressibility k = U^ 1 and a superfluid 
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FIG. 1: Possible renormalisation flows. The superfluid (SF) 
is stable for a large initial value of the inverse compressibility 
ft" 1 (= U). A small initial ft -1 may lead to phase separation 
(PS), while intermediate values give rise to a metal, especially 
at small doping. 



density n s = n. We are interested in how the interaction 
terms in £ lnt affect these quantities. The system will 
remain superfluid as long as both of these quantities are 
finite. 

Setting a = in the action above corresponds to the 
neutral superfluid. The resulting action is not quadratic, 
but neither of the remaining interaction terms modify 
the superfluid density (this must be equal to the total 
density for a Galilean invariant system). We therefore 
drop these terms, keeping only interaction terms that 
contain the gauge field, a. If the discarded terms have an 
effect on the compressibility then this can be absorbed 
into U. 

Any effect on the superfluid response must therefore 
come from the terms involving the gauge field. Simple 
perturbation theory is not suitable for discussing the in- 
stabilities at small n s and small U, since strong renor- 
malisation_pf these quantities renders it invalid. Feigcl- 
man et alB used a self-consistent method to discuss the 
reduction of the superfluid response in a system with 
long-ranged repulsion (so the compressibility is not renor- 
malised) . 

We use a different method, which allows us to treat 
the corrections to both the gauge field propagator and 
the bosonic propagator on an equal basis. We integrate 
out the highest frequency fields, which gives an effective 
theory for the remaining fields, with renormalised super- 
fluid fraction and compressibility. This process is then 
repeated until only the very smallest frequencies are left. 

As the bandwidth is reduced, three scenarios are pos- 
sible; these are shown schematically in figure [|. If both 
n s and n remain finite after all the fields have been in- 
tegrated then the system will remain superfluid. If k 
diverges at finite n s then there will be an instability to 
phase separation, but if n s vanishes at finite k then there 
will be an instability to a homogeneous non-superfluid, 
which we call a metal in the rest of this paper. 

To proceed, we generalise our model to allow for a 
superfluid density that is not equal to the mean boson 
density and a compressibility that is not equal to U . 



We write C x — C s x + C x , where 

+ |^|V^| 2 + \at- {Dl_ x ,y l £ (6) 

and 

4 = - Px (a x -\7e x ) + ^-(a^) 2 p x (7) 
m Am 

with ip%_ x ,)- x = (n s /m)S(x - x') + D~ l . Note that 
the phase stiffness of the bosons and the Meissner re- 
sponse of the gauge field are characterised by the same 
n s . This is required by the underlying U(l) gauge sym- 
metry. In the perturbational calculation the correction 
to the phase stiffness comes from the first term in C\ and 
the correction to the Meissner response comes from the 
second term. However, both corrections are given by the 
same polarisation insertion (the first bubble in figure^); 
this ensures that they are equal, as required. 

We start the renormalisation process with k _1 = U 
and n s = n. At each stage, we integrate out all boson ex- 
citations with frequencies whose magnitudes lie between 
the maximum u n and u n — 6Lu. We then average over fast 
fluctuations in the gauge field, as done in Refs. As 
each frequency shell is integrated out, we correct the val- 
ues of n s and k and use these values, denoted as n a (u> n ) 
and k(w„), for the next shell. 

The one-loop corrections to n s and n are given by: 

^ = -(<MII« , <5(«- 1 ) = -(5u>)K$ (8) 
n s 

where 
and 
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These depend on the gauge field propagator D (about 
which we have so far assumed nothing) and the propaga- 
tor for phonons in the Bose liquid: 

C LOn .q = — (PLU n ,qP—LU n , — q) 



w 2 +n g 2(«=i + _4_) 

These corrections correspond to the diagrams in fig- 
ure |^. To reiterate, the renormalisation scheme is im- 
plemented by using the renormalised values of n s and k 
when evaluating C Un and D UJn at each stage of the pro- 
cedure. Perturbation theory corresponds to performing 
the frequency integrals keeping n s and k constant. Our 
renormalisation scheme resembles a frequency dependent 
mean-field theory in that these parameters are constant 
with respect to momentum, but dependent on frequency. 
We see that Sn s is negative and 6k is positive: the per- 
turbations reduce the superfluid correlations and increase 
the density fluctuations. 
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FIG. 2: Polarisation insertions used in one- loop perturba- 
tional calculations. Dotted line is a gauge field propagator, 
solid line is a propagator for phonons in the Bose liquid. 
The processes affect the phase stiffness and Meissner response 
(left) and the compressibility (right). 

A. Overdamped propagator 

We now carry out this renormalisation scheme using 
the gauge field propagator given in equation || and in- 
vestigate the resulting phase diagram. Let us now iden- 
tify four independent dimensionless parameters for the 
system. We choose Aq = (ni/n)(m/m,{), U r — mil, 
m r = rrif/m and Q r = Q/E-p. Ao is the penetration 
depth of the Bose superfluid at n s — n in units of the 
fermion separation (i.e. lattice spacing at nt — 1). With 
the other dimensionless parameters fixed, varying Ao cor- 
responds to changing the boson density n: Ao oc 1/n 1 / 2 . 
U r measures the relative size of the repulsive and kinetic 
energies of the bosons. U r should be at least 0(1) if we 
want to model the hard-core repulsion in the slave-boson 
model. The parameter m r determines the coupling be- 
tween the bosons and the gauge field. The limit m r — > 
at constant Ao corresponds to the vanishing of the per- 
turbative corrections to C s . We expect our results to be 
reasonable for small m r . 

For convenience, we use a rescaled frequency, w — 

1/0 

wnjf/rif, and rescaled momentum u = pn t . We also 
define dimensionless propagators D WiU — nfD Untq /m{, 
and C w>u = ri{C Un , q /m%. 
The result is that: 

and 

(Mn u) = (*-k r ^ ^ (13) 

it J 2ir 

As mentioned above, we can see from these expressions 
that m r can be regarded as a measure of the coupling 
strength: the effect of the gauge field on the bosons is 
small if m r is small. Our choice of renormalised variables 
means that D s is independent of m r and U r . 

Performing the scaling process numerically (at Q / Ef = 
1) gives a phase diagram as shown in figure |3[ Our expan- 
sion (Q) about the superfluid state in means that these 
results are most reliable in the superfluid phase, up to 
the phase boundaries, which represent the onset of the 
instabilities of the superfluid state. The metal/phase- 
separation boundary in figure [3] is a rough guide only. 




FIG. 3: Phase diagram for m r = 1. SF: superfluid; PS: phase 
separation; Metal: vanishing superfluid response but finite 
compressibility. 
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FIG. 4: Position of critical point at varying m r . ac _1 and 
n s vanish together at the critical point U = Us, Ao = A3. 
Metallic region appears for m r > 0.29. 

We find that, at small U r , there is an instability to 
phase separation. This is because the correction to kT x , 
and hence U, from (||) is not directly dependent on U. On 
the other hand, at large U r , we find that the superfluid is 
stable. This is expected because IlW is small compared 
to unity and K^> is small compared to U, so that the 
corrections that we calculate are small. 

For intermediate values of U r , we find that there may 
be a metallic state: this is most likely at small doping 
(i.e., a large bare penetration depth Ao). However, be- 
cause is a fractional correction and not an absolute 
one, we note that it is not guaranteed that reducing the 
doping will result in a metal even as n — > (Ao — > 00). 

The form of the phase diagram depends on m r (see 
figure ^) . There is always a Meissner effect in the limit 
m r — > (at constant Ao). For small m r , only the in- 
stability to phase separation occurs. A metallic region 
appears in the phase diagram at m r ~ 0.3. Increasing 
m r above 0.3 favours both instabilities at the expense of 
the superfluid. 

In terms of doping, the boson density decreases from 
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left to right in Fig. ||. We see that the maximal doping for 
which the metal exists is approximately 10% for m r = 1 
and 5% for m r = 5, provided that U r is not too small. 



B. Propagating gauge field 

We now compare the results for an overdamped gauge 
field with the results for a propagating (Maxwell) gauge 
field, as used by Feigelman et alB. The propagator is 
D Un , q — g 2 (uj 2 +c 2 q 2 )~ 1 . As before, we rescale the boson 
mass and the penetration depth, using the parameters in 
the gauge field propagator. The relevant dimensionlcss 
parameters in this case are U r = mil as before; \ 2 M = 
(mg 2 /n) which is the rescaled penetration depth, and 
a = g 2 /(mc 2 ): the coupling constant for the gauge field. 
Rescaling frequency and momentum via w — (w„/<? 2 ), 
u = cp/g 2 , the equations |l2| and [l3| still apply, the only 
changes being in the form of D, and in the replacement 
of m r by a and Ao by Am- In this case the momentum 
integrals in equations |^ and |l^ can be done analytically, 
yielding: 
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where 7 = n(uj n ) [ui 2 + {n s (u n )g 2 /m)\ /(4nmc 2 ), B = 
w 2 K 2 (w„)/(4nn s (u; rl )), and X = \/\ — 4B are all dimen- 
sionless and are defined only for ease of writing. 
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Phase diagram for a Maxwell gauge field with 
: 1. SF: superfluid, PS: phase separation. 



tend to dominate over the loss of superfluid response, un- 
less n s disappears at rather a high frequency. The phase 
diagram for moderate coupling strength, a = 1, is shown 
in figure [|. The instability to phase separation is domi- 
nant at this coupling, so there is no metallic phase. 
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Phase diagram for a Maxwell gauge field with 
3, showing phase separated (PS), superfluid (SF) 
and metallic (Met) regions. The inset shows the narrow phase 
separated region between the metallic and superfluid regions. 



At larger values of the coupling, around a ~ 2, a metal- 
lic phase does appear; however, our one-loop perturba- 
tion theory is not reliable at these couplings. A phase 
diagram is shown in figure [| It is also clear from figure ^ 
that the metal requires very large values of the parameter 
mU to stabilise the system against phase separation. So, 
in contrast to the overdamped case, the parameters of the 
propagating field need to be pushed to rather unphysical 
values to see a metallic phase. 

Another consequence of the singularity in K^} is that 
k is very strongly renormalised near any metal-superfluid 
phase boundary. There is evidence that k _1 will always 
vanish before n s in a very small region near to the bound- 
ary: see figure ^. The phase-separated region between 
the metal and superfluid regions narrows with increasing 
Am- At large values of \m, an apparent superfluid- metal 
transition does appear, but the numerical calculations 
show a discontinuity in the compressibility at the tran- 
sition. This may be a symptom of an intervening phase 
separated region, which is too narrow to resolved numer- 
ically. 

As mentioned previously, the results of the renormali- 
sation scheme are likely to be most valid in the superfluid 
phase, and on the boundaries of that phase. The Maxwell 
case is distinct from the overdamped case in that K L 



(i) 



diverges as n s 



0, which has a strong effects on be- 



haviour of the system very close to the phase boundary. 



As discussed by Feigelman et alB, the propagating 
gauge field combined with short-range repulsion tends 
to cause phase separation. This is especially apparent at 
small n s : the correction diverges at small n s and u. 
This means that the instability to phase separation will 



IV. CONCLUSION 

To summarise, we have analysed the effect of both 
overdamped and propagating gauge fields on the ground 
state of a Bose liquid. In the overdamped case, we find 
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two different instabilities of the superfluid phase. Specif- 
ically, the penetration depth for the gauge field may di- 
verge. This is interpreted as an instability towards a 
homogeneous metallic phase. Alternatively, we find that 
its compressibility may diverge. We interpret this to be 
phase separation. We cannot comment on any spatial 
structures in the new ground state because our treat- 
ment cannot identify any characteristic length scale for 
this instability. A more careful treatment of lattice ef- 
fects and the compressibility at different wavevectors is 
needed. 

If the overdamped gauge field is replaced by a propa- 
gating one, then the instability to phase separation dom- 
inates over the divergence of the penetration depth for all 
systems with local repulsion. This is in agreement with 
Feigelman et alB. 

Applying these results to the slave bosons that describe 
the charge dynamics of the t-J model, both instabilities 
are found at dopings up to around 10%. For U r > 1, 
a metallic phase is favoured at these densities. (A large 
U r is natural for the hard-core slave bosons.) Although 
precise numerical values may be affected by the various 
approximations in our treatment, this Bose metal ap- 
pears to be in the regime in parameter space relevant to 
a non- Fermi-liquid description of the normal state of the 
cuprates. 
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APPENDIX A: SLAVE-BOSON DECOUPLING 
OF THE t-J MODEL 

In this section we review the derivation of a bosonic 
model of the form given in equation [l] from the fcr J model. 
Our notation follows that of Lee and Nagaosala. 

The t-J model is definedJn equation |[ We apply 
the 'slave-boson' decouplingaiJ to this problem and write 



Following Lee and NagaosaB we make the mean field 
decoupling Y,Mi*fjo) = Xoe la,j ■ We again neglect 
terms quartic in the boson operators and we arrive at 
the imaginary time Lagrangian: 

L r = ftr,r [ 9t ~ W + ia lr] fi<r,T + 

b lr [9r - Mb + iffl?, T ] bi, T - 

i 



+ c.c 



c.c 



(A2) 



where af is the real part of A;. The gauge invariance of 
the original electron operators under the transformation 
fi<7 —> fio-e j hi — > biC %Bi allows us to identify the ctij as 
the U(l) gauge field associated with this gauge symmetry. 

If we assume that the relevant fields are slowly vary- 
ing in space and make the continuum approximation, we 
obtain a Lagrangian density given by: 



r — 
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where the parameters mb and mf are effective boson and 
fermion masses respectively. They are related to the orig- 
inal t and J, and to xo- The real 2-vector field a x is de- 
fined by aij tT — (ri —rj) -ai^ ri+ ^ T . In making the con- 
tinuum approximation, we move from an original model 
with square symmetry to an isotropic model. 

The fermions and bosons interact only via the gauge 
field, so integrating out the fermions will give an effective 
propagator for the gauge field. The effective gauge field 
action is, to quadratic order: 



f irT b\ where / is a fermionic spin-± field and & is a g _J_ y a *f^ jjjO q y ^ q ) 



bosonic field. The bosons represent charged 'holes', and 
describe the charge dynamics of the system. Neglecting 
terms quartic in the bosonic fields (valid at small doping) 
we arrive at 



(Al) 



where the complex field A^ is a Lagrange multiplier en- 
forcing the constraint of no double occupancy. 



— l — Y a l( UJ n,q)^s{^n,q)a±(uj n ,q) (A4) 

where we have fixed the gauge such that V ■ a = 0; only 
the transverse component of a is non-zero and is denoted 
a± . The size of the system is L and the inverse tempera- 
ture is (5. The scalar part of the gauge field has a propa- 
gator determined by Ug(u) n , q) = (pf(u>, q)pf(u>, q)) where 
Pf(u),q) is the Fourier transform of the fermion density. 
At small energy and momentum Hg(cu n ,q) tends to a 
constant, as predicted by Thomas-Fermi screening: we 
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model this field by a local replulsion U. This is consis- 
tent with the introduction of the scalar field to imple- 
ment the constraint of exactly one fermion or one boson 
per site. The propagator of the vector part of a is deter- 
mined by U^(u> n ,q) = T^ v (J^*{uj n , q) J[(u> n , q)) where 

J f M (w„, q) is the Fourier transform of the fermion current in the rest of the particle-hole continuum. 



(B5) 



and T^ v = S^ v 



q 1 



selects the transverse part. 



APPENDIX B: POLARISATION OPERATOR OF 
AN ISOTROPIC FREE FERMI GAS 

In this section we discuss the overdamped gauge field 
propagator in a little more detail. We first state the 
form of the free fermion polarisation that we use and 
then discuss the effect of the cutoff 57. 

The current-current correlation function of the free 
fermions is given by 

(J?(u,q)J?(-u J ,-q)) = A x 

/ ^ E + §)"(p + §XB1) 

where G° is the single-particle Green's function for a free 
electron gas. 

The imaginary part of this propagator is non-zero 
only in the particle-hole continuum, where q{q — 2kp) < 
2mf \u>\ < q{q + 2kp). Within this region, we have: 



As mentioned in section II, we introduce an high- 
energy cutoff, = Ep, to the spectral function of the 
gauge field propagator. The effect of f2 on the fluxes 



4> 




FIG. 7: Plot of root mean square flux through an area R 2 . 
Flux is measured in units of the flux quantum; area in units 
of plaquette size. The fermion density is one per plaquette. 
Solid line is = 1000-Ef (increasing £7 has no effect at these 
areas. Dotted line is Q = -Bf. 



ImlTff (w,q) 
1 



3Trm{q 




fe=fe„ 



where k-p is the Fermi momentum and 



fcmin = max V fcp " 2m f |w|, 



m { \uj\ _ q 
2 



The real part of the polarisation is given by: 

, \2ra\u? + q 4 

Ren s {cj,q) = — — 

12irmiq z 

for 2m{ \u>\ < q(2kp — q), and by: 



ReU^(uj,q) = 



\2m\J 2 + q 4 
127rmfQ 2 



(B2) 



(B3) 



(B4) 



through various areas is shown in figure 0. The mean 
square flux through an area R 2 is calculated from: 



/ d 2 rdV((Vxo) rT (Vxa),.« T ) 



4 [% P dp 2 
Jo 2-k 



(B6) 



If no cutoff is used areas much smaller than one pla- 
quette still have large fluxes through them, which is not 
compatible with the underlying lattice in the t—J model; 
introducing the cutoff at the Fermi energy strongly re- 
duces fluxes through areas smaller than one plaquette, 
leaving the behaviour at larger lengthscales unchanged. 
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